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Abstract 

Monopole and photon contributions to abelian Wilson loops are calculated 

00. 

CN . using Monte-Carlo simulations of SU(2) QCD in the maximally abelian gauge. 



> 

in 



The string tension is well reproduced only by monopole contributions, whereas 
photons alone are responsible for the Coulomb coefficient of the abelian static 



o , 

potential. 

I. INTRODUCTION 

QCD confines color, but color confinement mechanism is not yet clarified. A promising 
picture is the dual Meissner effect due to condensation of some color magnetic quantity [|],@. 

x 

This picture is realized in the confinement phase of lattice compact QED [0-0]. Especially 
interesting are the following facts: 

1) A dual transformation can be done, leading us to an action describing a monopole 
Coulomb gas [[§,§f|8|]. Monopole condensation is shown to occur in the confinement phase 
from energy-entropy balance. 

2) The monopole contribution alone can reproduce the full value of the string tension |J. 
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In the case of QCD, there is a difficult problem. We have to find a color magnetic 
quantity in QCD. The 'tHooft idea of abelian projection of QCD JED] is very interesting in 
this respect. The abelian projection of QCD is to extract an abelian theory performing a 
partial gauge-fixing such that the maximal abelian torus group remains unbroken. After the 
abelian projection, SU(3) QCD can be regarded as a U(l) x £7(1) abelian gauge theory with 
magnetic monopoles and electric charges, 't Hooft conjectured that the condensation of the 
abelian monopoles is the confinement mechanism in QCD ||10|| . 

There are, however, infinite ways of extracting such an abelian theory out of SU(3) 
QCD. It seems important to find a good gauge in which the conjecture is seen clearly to 
be realized. A gauge called maximally abelian (MA) gauge has been shown to be very 



interesting ||i]]-[T3||. In the MA gauge, there are phenomena which may be called abelian 
dominance |12| , |l4[ . Especially the string tension in SU(2) QCD is well reproduced by Wilson 
loops composed of abelian link fields alone in the MA gauge. Moreover the monopole current 
fc M (s) can be defined similarly as in compact QED ||. It is shown in the MA gauge that 
the abelian monopoles are dense and dynamical in the confinement phase, whereas they are 



dilute and static in the deconfinement phase [|T^ . 

Recently we have derived an effective U(l) monopole action in the MA gauge and in 
SU(2) QCD |15|JT6[] . Entropy dominance over energy of the monopole loops, i.e., condensa- 
tion of the monopole loops seems to occur in the confinement phase if extended monopoles 
|17| are considered |15|JT6|] . After the abelian projection in the MA gauge, infrared behaviors 
of SU(2) QCD may be described by a compact-QED like U(l) theory with the running 
coupling constant instead of the bare one and with the monopole mass on a dual lattice. 

If the monopoles alone are responsible for the confinement mechanism, the string tension 
which is a key quantity of confinement must be explained by monopole contributions alone. 
This is realized in compact QED |J . The aim of this note is to show that the same thing 
happens also in SU{2) QCD by means of evaluating monopole and photon contributions to 



abelian Wilson loops. Preliminary results are reported by the present authors |T1| and other 
group 
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II. ABELIAN PROJECTION AND GAUGE INVARIANCE 



We adopt the usual SU(2) Wilson action. The maximally abelian gauge is given ]TT| by 
performing a local gauge transformation V(s) such that 



S,fl 

is maximized. Here 

U(s,ix) = V(s)U(s,fi)V- 1 (s + ^. (1) 

After the gauge fixing is over, there still remains a U(l) symmetry. We can extract an 
abelian link gauge variable from the SU(2) ones as follows; 

U(s,p) = A(s,p,)u(s,p,), (2) 

where u(s,fi) is diagonal and A(s, ft) has off-diagonal components. It is easy to show that 
the above fields u(s, p) and A(s,fi) behave under the residual U(l) transformation d(s) as 
an abelian gauge field and charged matters, respectively: 

u(s, fi) — > u'(s, ft) = d(s)u(s, fx)d){s + ft), (3) 
A(s, jl) -> A'(s, fi) = d(s)A(s, /t)rf f (s). (4) 

Let us repeat that a U(l) invariant quantity written in terms of the abelian link variables 
u(s, fi) after an abelian projection can be rewritten in a SU(2) invariant (but complicated ) 
form using the original link variables U (s, fi) 0]. The gauge function V(s) which maximizes 
R is a functional of U(s, ft). First study the transformation property of V(s) under any SU (2) 
transformation W(s), fixing the U(l) ambiguity of V(s) in some way. Considering that a 
gauge-fixed quantity does not transform, we see 

V(s) -> V w {s) =d(s)V(s)W- 1 (s) } (5) 

where d(s) (e U(l)) ensures the form invariance of V (s) and V(s) and is determined 
uniquely by V(s) and W(s). Using the definition (|I]) and the transformation property (0), 
we get 
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(U(s,jl)) w = d(s)U(s,ji)S(s + fi) (6) 

for any SU{2) transformation W(s). Hence all U(l) invariant quantities composed of U(s, jl) 
(and u(s,fij) are automatically SU{2) invariant after the abelian projection. 

As an example, consider a £7(l)-invariant lxl plaquette variable up(s,fi,v) composed 
of u(s,jl) alone. Using (|1|) and (0), we get 

u P (s, /i, v) = ^Ti(u(s, ft)u(s + ft, v)v)(s + z>, jl)u\s, v)) (7) 
= ±Tr(U'(s, p)U'(s + A, 0)U'\s + 0, jl)U'\s, z>)), (8) 

where U'(s,jl) is a modified link field defined by U'(s,fi) = B(s, jl)U(s, jl) and B(s, jl) = 
V^ ! (s)A ji (s, jl)V(s). Under an arbitrary SU(2) transformation W(s), we see from @, @ 
and © 

B(s,p,)^W{s)B(s,ji)W\s), (9) 
U'{s,jl) -> W(s)U'{s,fi)W*{s + p.). (10) 

Namely, U'(s,fi) transforms like the corresponding original link field U(s,fi). The SU(2) 
invariance of up(s, fi, v) is seen also from this property. Moreover, it is seen how different the 
abelian plaquette variable is from the full one. up(s,fi, v) in a different abelian projection 
corresponds to a differently modified full plaquette variable composed of U'(s,fi). 



III. MONOPOLE AND PHOTON CONTRIBUTIONS TO ABELIAN WILSON 

LOOPS 

We show an abelian Wilson loop operator written in terms of u(s,jl) alone after the 
abelian projection is given by a product of monopole and photon contributions. Here we 
take into account only a simple Wilson loop, say, of size I x J. Then such an abelian Wilson 
loop operator is expressed as 

W = exp{i y £ l J ll (s)0 lt {s)}, (11) 
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where J^s) is an external current taking ±1 along the Wilson loop and O^s) is an angle 
variable defined from u(s, fi) as follows: 



u(s,ji) 



(12) 



Since J^is) is conserved, it is rewritten for such a simple Wilson loop in terms of an an- 
tisymmetric variable M^s) as J u (s) = d'M^s), where d' is a backward derivative on a 
lattice. M pv (s) takes ±1 on a surface with the Wilson loop boundary. Although we can 
choose any surface of such a type, we adopt a minimal surface here. We get 

W = exp{- i -J2M llu (s)f IMU (s)}, (13) 

where f pu (s) = d^O^s) — dyO^s) and <9 M is a forward derivative on a lattice. The gauge 
plaquette variable can be decomposed into /^(s) = f pu (s)+27m pu (s) where /^(s) G [— 7r, 7r] 
corresponds to a field strength and n pu (s) is an integer- valued plaquette variable denoting 
the Dirac string Q Since M pu (s) and n pu (s) are integers, the latter does not contribute to 
Eq. (|13|) . Hence f pu (s) in Eq. (|j~3|) is replaced by /^(s). Using a decomposition rule 

= -E^( S - s'M^OVi//^ - ctM a/1 )( s ') 

+ - /3po- <9« <9q' M pcr ( s' ) ] , 

we get 

W =W 1 -W 2 (14) 

W 2 = exp{2vri kp(s) D (s ~ s')-e aPp(T d a M pa (s')}, 

where a monopole current k p (s) is defined as k^s) = (l/47r) e fWl p J d a fi3 1 (s) following 
DeGrand-Toussaint ||. D(s) is the lattice Coulomb propagator. Since /^(s) corresponds 



1 This condition /^^(s) G [— 7r, 7r] is applicable only in the case of defining a smallest l 3 monopole. 
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to the field strength of the photon field, ^(WSj) is the photon (the monopole) contribution 
to the Wilson loop. To study the features of both contributions, we evaluate the expectation 
values (Wi) and (W2) separately and compare them with those of W. 



IV. SIMULATIONS AND EXTENDED MONOPOLES 

The Monte- Carlo simulations were done on 24 4 lattice from (3 = 2.4 to /3 = 2.7. All 
measurements were done every 30 sweeps after a thermalization of 1500 sweeps. We took 50 
configurations totally for measurements. The gauge-fixing criterion is the same as done in 



Ref. |19fl . Using gauge- fixed configurations, we evaluated monopole currents and obtained 
the ensemble of monopole currents. 

As shown in the previous note [|l^,|l^], type-2 extended monopole loops with b > b c ~ 



5.2 x 10 3 (A^) 1 seem to condense, where b = na(/3) for n 3 extended monopoles and a(/3) is 
the lattice constant. So we measure 2 3 (3 3 ) extended monopoles with b = 2a(P)(b = 3a(/?)) 



of the type-2 |L7| as well as usual smallest ones. The extended monopole of the type-2 
is defined on an extended cube as the sum of the smallest ones included in the cube as 
shown in Fig. [l]. Note that the definition of the type-2 extended monopoles corresponds 
to making a block spin transformation of the monopole currents with the scale factor n. 
Hence the effective lattice volume is reduced. We call the effective lattice as a renormalized 
lattice. We evaluate the averages of W using abelian link variables (called abelian) on the 
original lattice, of W\ (photon part), and W2 (monopole part) on each renormalized lattice, 
separately. 

V. RESULTS 

The results are shown in the following: 

1. The monopole contributions to Wilson loops are obtained with relatively small errors. 
Surprisingly enough, the Creutz ratios %(/, J) of the monopole contributions having 
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small errors are almost independent of the loop size for example as shown in Fig. |2|. 
This means that the monopole contributions are composed only of an area, a perimeter, 
and a constant terms without a Coulomb term. 



2. Assuming the static potential is given by a linear + Coulomb + constant terms, we 
try to determine the potential using the least square fit. There are various ways, but 
we adopt a method similar to that pO using the Creutz ratios. The assumption of 



the form of the static potential leads us to the Creutz ratio 

x(i, J) = xo- xi(^i) + + X2( J(J _ 1) 1 J(J _ 1) ), (15) 

where xo is the string tension and xi corresponds to the Coulomb coefficient of the 
static potential. Using the fitted values of Xo an d Xi, we can reproduce each static 
potential. 

We plot their data in Fig. | (at (3 = 2.5) and in Fig. | (at (3 = 2.6). We find the 
monopole contributions are responsible for the linear-rising behaviors. When the 2 3 
extended monopoles are used, we obtain almost the same results. The photon part 
contributes only to the short-ranged region. There seems to exist a small discrepancy 
between the abelian and the monopole + photon parts for R/a = 12, but finite-size 
effects are expected there. Moreover, the assumption of the form of the static potential 
looks too simple. Similar data are obtained for (3 = 2.4 and 2.7. 

3. This is seen more clearly from the data of the string tensions which are determined 
from the static potentials. They are shown in Fig. [5[ The full and the abelian string 
tensions have large errors but they are seen to be well reproduced by the monopoles 
alone for /3 < 2.7 and the photon part has almost vanishing string tensions. 

4. The string tensions from the monopoles and the photons of various sizes are plotted 
in Fig. It is interesting to see that they are almost independent of the extendedness 
contrary to our preliminary data on a smaller lattice |?T] , although the monopole 
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actions determined in |15| , |T6l| depend on the extendedness. Note that the extended 



monopoles of the type-2 are composed of the sum of the smallest monopoles. Hence 
even when the extended monopoles alone look to condense, the smallest monopoles 
also have some information of the condensation. 

5. We have derived also Coulomb coefficients from the static potentials as shown in Fig. 
[7]. The monopole part has almost vanishing Coulomb coefficients which is in agreement 
with the constant behaviors of the Creutz rations of the monopole part as shown above. 
The l 3 photon part has large coefficients and they reproduce well the coefficients of 
the abelian static potentials. 

6. We also measured the Coulomb coefficients of the photon parts on each renormalized 
lattice. The values depend on the extendedness. But the following fact is interesting. 
The photon parts are evaluated on a renormalized lattice with b = na(f3). Hence they 
have different values of b for different extendedness. The Coulomb coefficients of the 
photon parts are well reproduced by the SU(2) running coupling constants g(b) with 
b = na(j3), i.e., — g(b) 2 /16n, where 

The scale parameter A determined is A ~ 48A^ which is quite near to the value 



A ~ 42Ai fixed from the monopole action [15 



In conclusion, our analyses done here strongly suggest that abelian monopoles are respon- 
sible for confinement in SU(2) QCD and condensation of the monopoles is the confinement 
mechanism if the abelian projection is done in the MA gauge. To extend our method to 
a finite-temperature system and also to SU(3) with or without dynamical quarks is very 
interesting. These studies are also in progress. 

We wish to acknowledge Yoshimi Matsubara for useful discussions. This work is finan- 
cially supported by JSPS Grant-in Aid for Scientific Research (c) (No. 04640289). 
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FIGURES 

FIG. 1. An extended cube on which an 2 3 extended monopole is defined as the sum of eight 
(=2 3 ) smallest monopoles. 

FIG. 2. Creutz ratios x(7, J) from abelian and 2 3 monopole Wilson loops at (5 = 2.6 versus 
I x J. The monopole Creutz ratio values are devided by 4, being adjusted to those in unit a([3). 

FIG. 3. Static potentials aV(R) versus R/a at (3 = 2.5. The values are shifted by a constant. 

FIG. 4. Static potentials aV{R) versus R/a at (3 = 2.6. The values are shifted by a constant. 

FIG. 5. String tensions at (3 = 2.4, 2.5, 2.6, and 2.7. 

FIG. 6. String tensions on each renormalized lattice. 

FIG. 7. Coulomb coefficients. 

FIG. 8. Coulomb coefficients on each renormalized lattice. 



11 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/94040 1 5 v 1 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/94040 1 5 v 1 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/94040 1 5 v 1 



This figure "fig2-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/94040 1 5 v 1 



0.1- 
0.08- 
0.06- 
0.04 
0.02 

0- 
-0.02- 
•0.04- 



1 1 1 1 1 r 

CREUTZ RATIOS beta=2. 6 



□ abelian 

# 2 3 monopole 



• • • • • • 



J I I L 



J L 



20 40 60 80 100 120 



I X J 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/94040 1 5 v 1 



This figure "fig2-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/94040 1 5 v 1 



0.8 



0.6- 



0.4- 



0.2- 



1 1 1 1 1 1 1 1 

beta= 

□ abelian 


1 


1 I I 




monopole 




$ 

s i 

f 

■ 




photon 




— 


monopole+photon 

■ 

■ 

a 






a 

1 

• 




■ 

■ 


• 


B 


■ 


■ _ 


B 

- • * 






a * * 


B 


m 

m 

m 






■ 

■ 


m 




i 


1 i i 



2 4 6 8 10 12 

R/a 



This figure "figl-4.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/94040 1 5 v 1 





1 1 1 1 1 1 1 1 1 | 1 1 

□ abelian beta=2.6 


0.5 


1 3 


mnnnnnlp — 

1 1 IVJI iuuuic 




— ▲ ■ 


nhoton 




• 1 


monopole+photon S m 






I 






I 






s 






s 






s 


0.25 




& 




3 


A ■ - 






■ 




• 

- B 


■ 

■ 

■ 






■ 






■ 

■ 




■ 


■ 

■ 





■ 


Ill 



2 4 6 8 10 12 

R/a 



This figure "figl-5.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/94040 1 5 v 1 



string tensions 



3000 



- [] 



2000- 



1000- 



0- 



o full SU(2) 
□ abelian 
■ l^monopole 
o 1 ^ photon 



2.4 



23 



2.6 



§ - 



2.7 



BETA 



3000- 



i 



2000- 



1000- 



string tensions 

■ 1 ^ monopole 
• 2 3 monopole 
a 3 3 monopole 



/ photon o 
/ photon o 
/ photon a 



4 



o- 



2.4 



2.5 



2.6 



2.7 



BETA 






I 1 1 
— + + 


1 i 1 i 

+ * - 


-0.1 




- 

— ■ 




Lambda=48(Lambda)L 




o abelian 


Coulomb coefficients 


-0.2 


■ 1 3 photon 






+ 1 3 monopole 
i , i 


i , i 



I I I I I I I 

2.4 2.5 2.6 2.7 

BETA 



~i 1 1 1 1 1 r 

o abelian * 2 3 photon 

■ 1 3 photon - 3 3 photon 




BETA 



